YALE-P8-00 
BUHEP-00-16 



The 2+1 Dimensional NJL Model at Finite Temperature 



Thomas Appelquist* 

Department of Physics, Yale University, New Haven, CT 06511. 

Myck Schwetz^ 

Department of Physics, Boston University, Boston, MA 02215. 

February 1, 2008 



Abstract 

We decribe properties of 2 + 1-dimensional Nambu-Jona-Lasinio (NJL) models at 
finite-temperature, beginning with the model with a discrete chiral symmetry. We 
then consider the model with a continuous U{1) x U{1) chiral symmetry, describing 
the restoration of the symmetry at finite temperature. In each case, we compute the 
free energy and comment on a recently proposed constraint based upon it. We conclude 
with a brief discussion of NJL models with larger chiral symmetries. 



1 Introduction 



The behavior of QCD at finite temperature and density has received renewed attention 
recently with the coming of the RHIC experimental program. The study of this phase 
structure is difficult for any strongly coupled quantum field theory, and it can be helpful to 
examine the problem in nontrivial but tractable models. In particlular, theories in less than 

four space-time dimensions can offer interesting and complex behavior as well as tractability, 
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often in the form of large-N expansions. In the case of three space-time dimensions, they can 
even be directly physical, describing various planar condensed matter systems. Models of 
interest for the study of finite temperature and density involve fermionic degrees of freedom 
and fall into two broad classes. One utilizes four-fermion interactions of the Nambu-Jona- 
Lasinio (NJL), or Gross-Neveu, form, and the other includes three dimensional gauge theories 
and closely related Thirring models. 

We restrict our attention in this paper to NJL models at finite temperature T. The 
features we discuss are interesting in their own right and will play an important role in 
understanding temperature-density phase diagrams for these models. All the analysis will 
make use of a expansion where Nf is the number of fermion species. The theories 

under consideration are renormalizable in the expansion unlike in the loop expansion 
. The expansion is reliable for all temperatures T except those in the vicinity of a certain 
critical temperature Tc. 

The question of interest is the spontaneous breaking of chiral symmetry. Although this 
is in some ways easier to handle in three space-time dimensions (3i5) than in four, there is 
an important subtlety in the case of continuous symmetry. The Coleman-Mermin- Wagner 
(CMW) theorem stipulates that the spontaneous breaking of a continuous symmetry can 
not happen in "iD at finite T [0]. This statement is related to the fact that it is impossible 
to write down a consistent theory of massless scalars in 2D. If a spontaneous breaking of 
continuous symmetry were to happen at finite T, then one would be faced with this problem 
at momentum scales below T, i.e. it would be impossible to construct an effective 2D theory 
of the Goldstone bosons' zero-modes. 

We first review the behavior of the NJL model with a discrete chiral symmetry, where 
the CMW problem does not arise. The broken symmetry at T = remains broken at finite T 
up to a critical value T^. We then turn to the 3Z) NJL model with a continuous U{1) x f/(l) 
chiral symmetry, exploring symmetry restoration at all T > 0. A critical temperature 
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continuous to exist, but it now marks a transition from the ordinary symmetric phase at high 
temperature to a low temperature Kosterhtz-Thouless phase that is also chirally symmetric. 
The broken phase exists only at T = 0. 

For each model, we also compute the thermodynamic free energy and enumerate the 
thermodynamic degrees of freedom. In a recent paper 0, the free energy F{T) was used as 
the basis for a proposed constraint on the behavior of asymptotically free theories. It was 
observed that for any theory governed by a fixed point in the ultraviolet (UV) or infrared 
(IR), the dimensionless quantity f(T) = — 27r[F(T) — F(0)]/,^(3)T^ approaches a finite value 
in the corresponding limit ( fuv = f(T oo) and fjR = f(T 0) ), and counts the 
(effectively massless) degrees of freedom if the fixed point is trivial. It was noted that 
fuv > fiR for all known asymptotically free theories in which the IR behavior is also free, 
or weakly interacting, allowing to be computed. It was conjectured that this is true 
for all asymptotically free theories. The theories considered in this paper are governed by 
nontrivial UV fixed points and apriori may or may not satisfy this condition. 

The outline of the paper is as follows: in Section 2, we review the NJL model with a 
discrete chiral symmetry; in Section 3, we treat the NJL model with a continuous U{1) x 
f/(l) chiral symmetry; we summarize and briefly describe NJL models with larger chiral 
symmetries in Section 4. 

2 Discrete Chiral Symmetry 

The NJL model with a discrete chiral symmetry for Nf copies of Dirac fermions is described 
by the following LagrangianQ: 
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^We adopt the notation of for 7-niatrices: 7'' ~ ® 
dimensfonal Minkowski space with (1, —1, —1) signature. 
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where ^' may be taken to be a 4A^j-component fermion field. The discrete symmetry is 
^ — > 75^, and Qq is the couphng. We analyze the model using the 1/Nf expansion. 

We work first to leading order and then note that for this discrete-symmetry model 
the qualitative behavior is not modified in higher orders. It is convenient to introduce an 
auxihary field a coupled to The finite temperature effective potential as a function of 
a may be computed to this order by integrating out the fermions: 

2 

Veff{<7)/Nf = + ^Trln{^<^-a), (2) 

with the second term evaluated at finite temperature. 

The only cutoff dependence is in the renormalization of go, and it may be carried out 
at zero temperature. Extremizing the zero-temperature effective potential {dVeff/da = 0) 
gives 

where A is an ultraviolet cutoff and where the extremal value ctq describes the fermion mass 
gap at zero temperature. As (?gA approaches the critical value 7r^/2 from above, Uo/A 0. 
Equivalently, the cutoff may be removed holding ctq fixed providing that g^A is tuned to 
7r^/2. With this prescription, it can be seen that the high energy behavior of the theory is 
described by a nontrivial but weak ultraviolet fixed point. That is, a finite effective four- 
fermion interaction g^{p)/2Nf is induced such that for p >> ctq, the dimensionless quantity 
g\p)-P^O{l). 

To leading order in the zero-temperature effective potential is now 

Vl,fW)IN, ^ - "^a^ (4) 



which has a minimum at cr = ctq and is convex for cr > cro/2. At finite temperature, the 
leading order effective potential takes the form: 

V.„{o)IN, = l-a^' - - r f'^lnil + eM-f^]) ■ (5) 



It may be extremized to give the temperature dependent gap equation 

(XT = - 2Tln{l + exp[-^]) = ao - 2Tj2^^ exp[-k^] . (6) 

This leads to the critical temperature T^. = cro/(2/n2), above which cry vanishes. The 
existence of a finite for a discrete symmetry in 2 + 1 dimensions is perfectly acceptable 
and this qualitative feature is not changed by the higher order terms in the 1/iV/ expansion. 

Before discussing higher order terms, we comment on the thermodynamic free energy of 
this model which is given by F{T) = Veff^ax)- Then, to leading order in the /(T) 

defined in the introduction will take the form 

J V ; ((3)T'^ V 3 2^6 A;^ V ) f[ \ ) , \ ) 

where ax is determined by the gap equation (P). It can be shown that this function increases 
monotonically from //r = at T = to fuv = 3A^/ at T = oo. The infrared value, fm = 0, 
is expected because a mass gap develops for for T < Tc. The ultraviolet value, fuv = 3iV/, 
counts the number of fermionic degrees of freedom (4 Nf) times the Fermi-Dirac factor 
(3/4). This is expected because the model is governed in the ultraviolet by a weak UV fixed 
point - there will be higher order corrections in the 1/iV/ expansion. Monotonicity may be 
established by noting that 

df{aT,T) ^ d^doT , dl ^ 
dT dar dT dT 

-^(4-f44).^^Kl.e..|-^,). (8) 

The first term is always non-negative since ax < ctq; the second term is always positive. 

Finally we describe the effect of the higher order terms in the expansion. Through 
next order, the effective potential is given by 

2 

V,fj{a)/Nj = ^ + iTTln{i<^-a) + ^Tr In D-\ (9) 
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where D^^{p^) is the inverse a propagator computed in a background a field at finite tem- 
perature T. At zero temperature, this propagator is given by 

Using the leading order gap equation (|]) and taking the limit A — > oo with ctq fixed, D^^{p^) 
takes the form 

Z^;i(p2) = ^ a-ao + \ tan-'^] . (11) 

For nonzero temperature D^^{p'^) is given in Ref. 0. Expressions, similar to Eq. ^may be 
written down for higher order terms. 

For T well above Tc, the 1/A^/ expansion for Veff{cr) may be seen to converge for large 
Nf. As T ^ Tc from above, however, T — T^. In this limit, higher order terms in 

Veff{(y, vr) become singular and trigger the breakdown of the 1/A^/ expansion. They may be 
described by an effective 2D, Landau-Ginzburg theory, consisting of the zero mode of the a 
field, relevant at scales below T. The terms in the Landau-Ginzburg Lagrangian, in addition 
to the common mass term are a kinetic term, and a Acr^/4! interaction term. Each 

may be computed to any order in l/Nf by integrating out the fermions. To leading order, 
each arises with coefficient Nf. 

Using this effective 2D theory and counting infrared powers, it may be seen that the 
effective dimensionless expansion parameter for V"e//(cr), for small a, is of the form 

1 r2 

(12) 



where = m^ + Ao"^, up to logarithmic corrections. Thus in the neighborhood of the origin, 
the 1/A^^ expansion breaks down for |T — T^/T ^ 1/A^/. In particular, it breaks down for 
the free energy and f{T). This is true also as T approaches Tc from below, as indicated by 
the absolute value sign 0. 

The 1/A^/ expansion remains convergent as long as \Tc — T\/T » l/Nf. In the high 
temperature limit, it correctly describes the effective potential and free energy in the sym- 
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metric phase, leading to small corrections to the leading order estimate fuy = SNf. In the 
low temperature limit, it describes the broken phase, and because of the mass gap, leads to 
fiR = to all orders. The inequality fuv > fm is satisfied. 



3 Continuous U{1) x U{1) Chiral Symmetry 

The NJL model with a continuous U{1) x f/(l) symmetry \E' — > explia'j^)'^ , \E' exp{if3)'$ 
is described by the Lagrangian 0: 

2 

£ = + ^[(M)' - (%5^)']. (13) 

The model may be analyzed to leading order in the l/Nj expansion by introducing auxil- 
liary fields a and vr and integrating out the fermions. The leading-order, finite-temperature 
effective potential takes the form 

2 2 

V,ff{a,7i)/Nj = + zTrln U^-a-zj.n) (14) 

As in the case of discrete symmetry, the renormalization may be performed by extremizing 
the effective potential {dVejf/da = and dVeff/dn = 0) at zero temperature. The vacuum 
expectation values may be rotated using the chiral symmetry so that tt has a vanishing 
expectation value. The zero temperature gap equation then takes the same form as in the 
discrete-symmetry model: 

l-^^W(l->,-'A), (15) 

where A is an ultraviolet cutoff. As before, Uo may be held fixed in the "continuum" limit 
A — i> oo {g^A —>■ vr^/2) describing the fermion mass gap at zero temperature. The chiral 
symmetry is broken and the vr field describes the associated Goldstone boson. The zero 
temperature effective potential to leading order in 1/A^/ now reads: 



which has degenerate minima at + tt^ = cTq and is convex as a function of two variables 
if + > al 

At finite temperature, the leading order effective potential has the same form as Eq. (^), 
with the replacement cr^ ^ + vr^: 

vr)/iV, = i- {a-Wf^ - ^ {a-W) - f ^ ln{l + exp[- f7^^^^]) . 

(17) 

Extremizing it suggests that, as in the case of discrete symmetry, the broken symmetry of 
the zero temperature theory remains broken at finite temperatures below a critical value 
Tc = o'o/2 ln2. But unlike the discrete case, this conclusion cannot be correct since the 
zero-mode of the associated Goldstone boson would describe an effective 2D theory with a 
spontaneously broken continuous symmetry - in contradiction with the Coleman- Mermin- 
Wagner theorem. 

We explore the resolution of this problem by first noting that as in the case of discrete 
symmetry, the 1/iV/ expansion is convergent as long as T is not near the transition temper- 
ature Tc ( \Tc — T\/T » 1/Nf ). The next order term in Veff{<J,n)/Nf, for example, may 
be written in the form 

In D;' + ^Tr In D;\ (18) 

where and are functions of T, T^, momentum, and field strength 0]. This term can 
be seen to be of order 1/A^/ for \Tc — T\/T » 1/Nf. But as in the case of discrete symmetry, 
the 1/Nf expansion breaks down due to infrared singularities when \T — Tc\/T ^ 1/Nf, with 
the singular terms describable by an effective 2D Landau- Ginzburg Lagrangian. 

For T » Tc, the theory is in the symmetric phase with a convergent 1/A^/ expansion. 
For T << Tc, even though the 1/Nf expansion is again convergent, the continuous symmetry 
model is, unlike the discrete symmetry model, not in the broken phase. The chiral symmetry 
remains unbroken, although in a way different than in the high temperature range. To see 
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this, it is convenient to use the following parametrization^ of the auxiliary fields for T < T^, ||^: 

0- + zvr = pe^^ . (19) 

The leading order potential, Eq. [1^, then depends only on p, and extremization leads to a 
non-zero VEV pr, equal to po = ctq at T = and vanishing like (Tc — T)^/^ as T ^ Tc- The 
fermion has mass pr, the fluctuations of p have mass 2pT and there is a massless scalar field 

e. 

This behavior implies symmetry breaking, however, only if 9 takes on some fixed and 
non-zero vacuum value. Now as is well known, this doesn't happen in 2D since quantum 
effects generate logarithmically infrared divergent fluctuations. In the present model at finite 
T , the same is true since the long distance behavior is effectively 2D. To explore this in 
detail, we first restrict attention to T sufficiently below ((T^ — T)/T » l/Nf) so that 
the expansion is a reliable tool . The realization of the symmetry may be studied by 

examining the theory in the infrared, at momentum scales well below pr, where the fermions 
and the p fluctuations may be integrated out. The only massless degree of freedom is 6, 
and if pr and T are of the same order, the non-zero Matsubara frequencies of 6 may also be 
integrated out. The effective low energy theory describing physics at momentum scales well 
below T and pt is then a 2D chiral Lagrangian p, ^: 




where the dots indicate higher derivative terms in 6. The dimensionless coefficient 1/t is 
given to leading order in the 1 / Nf expansion by 




^ This transformation of variables must be handled carefully. The jacobian of the transformation in the 
functional integral leads to a cubically divergent term in the effective lagrangian of the form (0) In p, which 
serves only to cancel other such terms arising in the quantum computation Without this cancellation, 
the jacobian term might be assumed pOj to change the cutoff dependence from that described in Eq. n3, 



eliminating the existence of the continuum limit as described there. 



9 



In the very low temperature limit (T << T^), physics at momentum scales of order T 
as well as below it may be described by an effective theory obtained by integrating out the 
fermions and the p fluctuations, but keeping all the Matsubara freqauncies of 6. This theory 
is described the effective chiral Lagrangian (|2QD, but in three dimensions rather than two. 

Since the underlying theory is Abelian, this effective Lagrangian describes a free massless 
scalar field 6. Chiral symmetry breaking depends on the behavior of correlation functions 
of physical, and therefore single valued functions of 6 {cos 6, sin 6, or equivalently e^* ^). 
This behavior is determined the coefficient t in the chiral Lagrangian, which is of order 
or smaller unless pr/T is small. (This only happens if T approaches Tc and we are 
avoiding that limit now to insure convergence of the 1/iV/ expansion). Analysis of this 
2D theory, including the role of vortex solutions reveals that for small t, it is in the 
Kosterlitz-Thouless phase where the contribution from vortices is negligible, and where, 
for example, the correlation function (e*^*-^^ g-«^(o)^ j^g^g characteristic power-law behavior 
(^Q^s{x) ^-ie{o)'^ ~ x^*/^'^ at distances x >> 1/T. The power-law fall-off, while corresponding 
to an infinite correlation length, still indicates an absence of long range order. There is no 
chiral symmetry breaking. 

The transition from broken to unbroken chiral symmetry is at T = 0. In the limit 
T — > the parameter t — >■ 0, and the range of relevance for the effective 2D theory and the 
(weakening) power-law fall off moves off to infinity. In the zero temperature 3D theory, the 
6 field develops a fixed VEV and describes a Goldstone boson. 

At finite T, the analysis of the effective Abelian 2D theory leads to the conclusion 
that when t becomes of order unity, the interactions become strong, vortex solutions play 
an important role, effectively renormalizing t, and a finite correlation length develops. This 
is the normal phase corresponding to in-tact chiral symmetry. For the U{1) model being 
discussed here, t can become of order unity only when p^/T ^ 0, that is, when T ^ T^. from 
below. To be more precise, since pt ~ 

(rc-T)i/2 in this limit, t becomes of order unity only 
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when (Tc — T)/T^ 1/^/- But we have aheady noted that this is the range where the 
expansion breaks down. The p-field becomes hght and must be included along with 6 in an 
effective 2D Landau- Ginzburg Lagrangian describing the transition at Tc, now interpreted 
to be the transition from the Kosterlitz-Thouless phase to the symmetric phase with finite 
correlation length. The infrared singularities associated with the 2D description mean that 
the 1/Nf expansion is not directly useful to describe this transition. 

Finally, we comment on the behavior of the thermodynamic free energy F{T) = Veffipr) 
for this continuous-symmetry model, and the proposed inequality constraint using the 
quantity /(T) = -27r[F(T) - F{0)]/C{S)T^. In the limit T ^ oo, the model is weakly 
coupled (as was the discrete-symmetry model) with the dynamics governed by a nontrivial 
UV fixed point of strength l/Nj. To leading order in 1/Nf, f{T) takes the same form as in 
the discrete-symmetry model: 

Thus to this order, fuv = 3A^/. There will be 0{1) corrections to this result due to the 
weakly interacting UV fixed point. 

In the low temperature limit (T ^ Tc), the discrete- and continuous- symmetry mod- 
els behave differently. In both cases, the expansion for the free energy is convergent, 
describing in the discrete case a mass gap and leading to fm = to all orders. In the con- 
tinuous case, there is one massless degree of freedom described by the 6 field. At momentum 
scales on the order of T which determine the free energy, its behaviour is fully 3D, governed 
by an effective 3D chiral Lagrangian of the form of Eq. Since this 3D Lagrangian is 
non-interacting, fjR may be computed exactly to give fiR = 1 corresponding to the single 
massless degree of freedom. This computation of fm is, in effect, an evaluation of the next 
to leading term in the Veff given by Eq. |18|, for T -C Tc. For this model as for the discrete 
symmetry model, fm < fuv- 
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4 Summary and Discussion 



We have described the finite temperature phase structure of two 3D NJL models analyzed in 
a expansion where Nf is the (large) number of four-component fermions. In the more 

familiar case of a discrete chiral symmetry, broken at T = 0, there is a transition from broken 
symmetry to unbroken symmetry at a non-zero temperature Tc. The low temperature phase 
is characterized by a mass gap; since there are no massless degrees of freedom, the quantity 
fiR defined in the Section 1 takes the value fiR = 0. The 1/Nf expansion breaks down as 
T ^ Tc due to infrared singularities describable by an effective 2D Landau- Ginzburg theory. 
At high temperatures, the 1/Nf expansion is again convergent for large Nf and the theory 
is governed by a weak 0(1/ Nf) UV fixed point. One finds fuv = 3A^/ up to corrections of 
0{l/Nf). 

In the case of a continuous f/(l) xf/(l) symmetry broken at T = 0, there is no spontaneous 
breaking of the chiral symmetry at T 7^ 0, although a phase transition still exists at a 
non-zero temperature Tc. The normal high-temperature symmetric phase changes at to 
a low-temperature phase with its low momentum components {p <^ T) described by an 
effective 2D theory in the Kosterlitz-Thouless (K-T) phase and with a power-law behaviour 
of correlation functions for x ^ 1/T. In both finite-T phases, the continuous chiral symmetry 
is unbroken. The transition from the K-T phase to the broken phase is at T = 0. There is a 
single massless degree of freedom for T < which becomes a conventional Goldstone boson 
in the limit T ^ 0. As in the case of discrete symmetry, the 1/Nf expansion breaks down at 
the phase transition T due to 2D infrared singularities, but is well behaved away from 
the transition for large Nf. We have nothing to say in this paper about the phase structure 
of NJL models at small Nf Q- 

The quantity fiR defined in the Section 1 takes the value fiR = 1, refiecting the presence 
of a single massless degree of freedom below Tc. The UV fixed point leads to the value for 
fuv = SNf up to corrections of 0{1/Nf) as in the discrete case. 



12 



For both models, the exploration of the transition at T = requires methods that are 
non-perturbative in 1/Nf. For the discrete case, e-expansion methods (4 — e dimensions) 
may be brought to bear although convergence is problematic since e = 2 . In the continuous 
case, the conventional e expansion is not useful since the model is in the broken phase for 
T < Tc for any e < 2. In both cases, the breakdown of the 1/iV/ expansion is associated with 
finite-temperature infrared divergences characteristic of a theory in less than four space-time 
dimensions, and has no direct counterpart in four dimensional theories. 

It is interesting to generalize our discussion to NJL models with larger continuous chiral 
symmetries. Suppose, for example, that the continuous U{1) ® U{1) symmetry of ( pTSD is 
extended to a U{n)®U{n) symmetry which at zero temperature breaks to the diagonal U (n) 
producing Goldstone modes. For the case << Nj, the expansion may be used to 
analyze the symmetry breaking pattern at zero and finite temperatures. An order by order 
analysis concludes that for T > Tc ^ po (where po is a zero temperature vev of an auxilliary 
field) the model will be in the symmetric phase. As in the abelian case the expansion 
breaks down when T ~ Tc due to 2D infrared singuliarities. For T < Tc, the fermions are 
massive, as in the abelian model, and may be integrated out to describe lower energy physics. 
At momentum scales below T, an effective 2D chiral Lagrangian emerges, which naturally 
breaks into two parts. One is the free abelian model we considered previously and the other 
is an SU{n) x SU{n) chiral Lagrangian. The latter is interacting and asymptotically free 
i, 0,0. 

The abelian model will be in the KT phase as long as T is not close to Tc, ensuring 
that the coupling strength is small - of order l/Nj. At intermediate momentum scales, 
p < T, the coupling strength of the non-abelian part is also small - of order l/Nj. The 
/9-function for this coupling t is (3}: —ht'^ where the constant h is positive and of order 
n^/(47r) for n> 1. Thus the effective i coupling runs, becoming of order unity at ultra-low 
scales ~ T exp(— 27r ■ Nf/in? ■ pr/T). As in the case of the U{1) ® U{1) model, unit 
strength is expected to disorder the system leading to a finite correlation length, now of 
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order But since this happens for any finite T < T^, the non-abehan theory is in 

the ordinary symmetric phase, not the KT phaseQ For the abehan part only, there is a 
KT transition as T — Tc. The value of the fm is because although the light scalars of 
the non-abelian part of the effective theory are massive, their mass msu{N) ~ f^m vanishes 
exponentially as T — > 0. Thus they are effectively massless in this limit. Clearly, for generic 
<^ Nf the inequality fiR < fuv will always be satisfied since fuv = 3 A^/. 

If the U{n)®U{n) model is arranged so that becomes of order Nf, the expansion 
breaks down at all T due to the large number of scalar degrees of freedom {0{n'^)) that are 
formed. It is beyond the scope of this paper to analyze the model when ~ Nf, but we note 
that it may well be accessible to a combined '^/Nf, l/n^ expansion. In the limit T << px, 
the model is described by an effective, 3D chiral Lagrangian at the scales of T, and the 0{N) 



models of this sort have been analyzed using the expansion |l5 . 
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